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The complex structure of a surface generated by the two-dimensional dynamical triangulation(DT) is deter- 
mined by measuring the resistivity of the surface. It is found that surfaces coupled to matter fields have well-defined 
complex structures for cases when the matter central charges are less than or equal to one, while they become 
unstable beyond c = 1. A natural conjecture that fine planar random network of resistors behave as a continuous 
sheet of constant resistivity is justified numerically for c < 1. 



1. Introduction 

The discretized approach of 2D quantum grav- 
ity, implemented by the matrix model, exhibits, 
in a continuum limit the same behavior as the 
continuous approach, given by the Liouville field 
theory. Thus it is certain that the two theories 
describe the same physics. However, we feel that 
the mutual relationship is not yet fully under- 
stood. While the relation between DT and matrix 
model is transparent, its relation to the Liouville 
field theory is not so evident. Intuitively, we ex- 
pect that the simplicial manifolds generated by 
DT will tend to continuous manifolds in the large 
number limit of simplexes. However, such sur- 
faces are known to be fractal [Q, hence the way 
the lattice regularization simulates the continu- 
ous theory is not at all clear. It is therefore indis- 
pensable to make sure that an universal contin- 
uum limit exists. In order to check this numer- 
ically, we propose to measure how the network 
conducts currents, and compare it with the corre- 
sponding continuous medium. An invariance un- 
der local scale transformation plays an essential 
role in the algorithm, and because of this prop- 
erty we can extract information about the com- 
plex structure which is independent of the con- 
formal mode parametrised by moduli. 

Despite many theoretical developments the 
precise nature of the transition at c = 1 is still 
not fully understood. Serious numerical efforts 
have failed so far to convincingly establish the 
transition, at least by measuring the string sus- 
ceptibility. It should be much easier to observe 
this transition in local scale invariant quantities, 



rather than in the original metric which is com- 
paratively complicated reflecting the fractal na- 
ture of the surface. Here, we show such a method 
by which we can define the complex structure of 
a DT surface and separate the conformal mode. 

2. Determination of the complex structure 

The basic idea of our method originates from 
the analog model || of the dual amplitude in 
analogy with the electrical circuit || . Let us first 
consider a two dimensional conducting medium 
with conductivity tensor density a^ u . Ohm's law 
reads 

f = a^E v , (1) 

where E„(= —d v V) is an electric field and is 
a current density. The Joule heat Q can be esti- 
mated by a potential distribution V on the sur- 
face; Q = / tfxa^d^VdyV. Requiring SQ = 0, 
it leads to an equation for V: d^a^d u V = 0, 
which is just the equation of continuity. For a 
homogeneous isotropic medium the conductivity 
tensor density is given by a^ iV = \^fgg^ v with r 
being the resistivity constant. Therefore, Ohm's 
law(|l|) is invariant under the local scale trans- 
formation g^ v i — ► g Al „e~ <T . In the following we 
first give an algorithm to measure the resistivity 
from voltage distribution. In the case of spheri- 
cal topology we can regard the surface essentially 
as an infinite flat sheet. Then the potential at a 
point specified by the complex coordinate, z, with 
a source of current / placed at Zk , and a sink of 



the current at z; is given by 

V(z) = -^ln|^— ^| + Const., (2) 

Z7T Z — Z[ 

In order to avoid the ambiguity arising from 
Const, in eq.(|^) we measure the potential drop 
between coordinates Zi and Zj . We define the po- 
tential drop between points i and j in the presence 
of a source and a sink of current at points k and 
Zby 

V# )W = V(z t )-V( Zj ) 
Ir 

= --^^m,^j;zk,zi]\, (3) 

where [zi, Zj-,Zk, zi] = z ^~-z z J -z' * s known as the 
anharmonic ratio. It has an important property: 
i.e. SL(2,C) invariance. Then it allows us to 
fix three coordinates among the {z^ Zj, zj,, z{\ to 
any desired values without changing the potential 
drop eq.([s]). For example, we can fix three coor- 
dinates as Zj = 1, Zfc = 0, and zi — oo. Then the 
potential drop eq.(||) for I = 1 is written as 

^f (0 =-^ ln NI. (4) 

V^=-^ln\l- Zi \. (5) 

All the other possible permutations of {z} only 
give the potential drops obtained by linear com- 
binations of eq.([|) and eq.(^|). Since unknown 
variables are {zj, r}, we need the potential drop at 
an additional point z m . Then we get three more 
relations in addition to eq.(^) and eq.(||) which 
altogether are sufficient to specify five unknowns: 
{zi,z m ,r}. The three extra equations are given 



by 

Vif l) =-^\z m \, (6) 

V^ (l) =-^ln\l-z m \, (7) 

Vt m +V<§®=~ln\z i -Z m \. (8) 



Now we apply this to a random surface generated 
by DT. Since the local scale invariance is not fully 
realized on this surface, the resistivity may have 



certain fluctuations. In fact, the value varies from 
one sample to another, and depending on the po- 
sition of electrodes for the measurements. Such 
ambiguities should, however, vanish in the con- 
tinuum limit, if it exists at all. The dual graph of 
a surface consisting of N triangles is regarded as 
a trivalent network, where we fix the resistance of 
the link connecting two neighbouring vertices to 
be 1. What we want to do is to examine whether 
such a network behaves like a continuous medium 
in the limit of N — > oo. For the measurement we 
pick five vertices in the dual graph and perform 
the method as explained above. The practical 
method we employ for the determination of po- 
tential drops is as follows; we pick two vertices at 
Vi n and v ou t for the source and the sink of cur- 
rent with unit intensity. By writing the poten- 
tial of the vertex v as V(v), current conservation 

reads V(v) = |{5^=i,2,3 V ( v i) + ~ 5 v,v out }, 
where v± , «2 and are the three neighboring ver- 
tices of v. We solved this set of equations by the 
successive overrelaxation method. 



3. Analog model 

Our method can be reconsidered from the 
string-theoretical point of view. Whether string 
theory can be described by local fields or not is 
one of the fundamental questions. The authors 
of P| showed that a sum of planar fishnet Feyn- 
man diagrams with N external particles is indeed 
approximated by a N particle Veneziano ampli- 
tude. The basic assumption of this model is that 
the network corresponding to a fine planar fishnet 
Feynman diagram can be regarded as a homoge- 
neous isotropic conducting sheet with a constant 
resistivity. Once this is accepted it is straightfor- 
ward to derive dual amplitudes by evaluating the 
heat generated by the sheet. The purpose of our 
article is to justify this basic assumption numer- 
ically. In order to clarify the correspondence to 
the network of resistors, we write the amplitude 
of fishnet Feynman diagrams as 

J II daip(ai) J J] {Pane-* (9) 

planar Zrpropagator i:vertex 
diagrams 



where Q = 2 a t > ^ an< ^ (*'•?') indicate a link 

and its two ends(vertices) of the network respec- 

_D j 2 

tively, and p(a t ) = 2(27r)D/2 a ; 2 e"2 Qim . Here 

we have introduced a Feynman parameter a; [j] 
for each propagator I in the momentum integra- 
tion, 

f dDp 1 Jp(*-v) 
J (2Tr) D p 2 + m 2 

r°° 1 _D i a (*~«) 2 

The basic conjecture is that the fishnet 
diagrams behaves as a continuum homoge- 
neous conducting sheet after taking the average 

E planar JILpr opagator 

daip{ai). 

diagrams 

4. Numerical results and discussions 

We fix the topology to be two sphere^ 2 ) and 
the number of triangles to be 4,000, 8,000 and 
16, 000 (prohibiting tadpole and self-energy dia- 
grams). For matter fields coupled to gravity we 
put n scalar fields with n — 0, 1, 2 and 3 and/or 
Ising spins on each triangle. 

4.1. Resistivity distributions 

We pick-up five vertices randomly as five elec- 
trodes to obtain the resistivity constant as indi- 
cated in Sec. 2, and repeat this procedure 50 
times for each configuration, for about 100 inde- 
pendent configurations as ensembles. 

Let us first take a look at the results for pure 
gravity (Fig.0). We find that peaks get narrower 
as the size grows. The value r ss 2.6 should be 
compared to V3 of the flat network (i.e. a hon- 
eycomb lattice). An increase in r is understood 
to be a reflection of the fractal nature of the sur- 
face. As the surfaces get larger, finite size effects 
due to the discreteness diminish and eventually 
the peak will grow indefinitely. This is what we 
have expected as the continuum limit of a net- 
work of resistors. We find similar tendency in the 
one scalar case(c = l)(Fig^). Also, in Fig^, we 

1 There are analogies between Feynman diagram and net- 
work of resistors: (a; <-> resistance), (xj <-> voltage at i), 
(pi <-» current through 1) and (Q <-> heat generation). 



show data to justify the assumptions of the ana- 
log model for the case of pure gravity. Each resis- 
tance is distributed randomly with weight e~° l . 
These distributions indicate the same behavior 
with previous cases c = 0. In Fig.|],[| we plot the 
distributions of r for the c = 2, 3 cases. Here no 
sharpening of peaks with increasing sizes is seen 
unlike the previous cases, which suggests that the 
surface does not approach a smooth continuum 
limit. This is considered as the reflection of the 
tachyonic mode instability in the continuous the- 
ory. The same tendency is also supported in the 
case of 3 critical Ising spins. 

In conclusion we have observed the so called 
c = 1 barrier for two dimensional quantum grav- 
ity coupled to matter as the transition from a 
well-defined complex structure to an ill-defined 
one. 




Figure 1. Resistivities for the case of pure grav- 
ity. The size of (a), (6) and (c) are 4000, 8000 
and 16000 triangles respectively. 
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Figure 2. Resistivities for the case of 1 scalar. 

4.2. Separating the conformal modes 

Next, we discuss how to extract information 
on the conformal mode. Once we have found the 



Figure 3. Resistivities for the case of pure grav- 
ity when each resistance is randomly distributed. 



Figure 5. Resistivities for the case of 3 scalars. 
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Figure 4. Resistivities for the case of 2 scalars. 

value of resistivity, we can assign a complex co- 
ordinate to each vertex of the dual graph from 
eq.(||) and eq.(|J). In Fig. 6 we plot the distribu- 
tion of coordinates for all the vertices of a sur- 
face for the case of pure gravity with 8000 trian- 
gles. The point density around a point z should 
be proportional to y/g(z), because each vertex 
is supposed to carry the same space-time vol- 
ume. In the conformal gauge which we have been 
employing, we have y/g(z) = exp0(z), where 
4>{z) is the conformal mode. This argument may 
be over-simplified in the sense that we did not 
pay attention to the quantum fluctuation of 4>{z) 
and the renormalization of the composite opera- 
tor exp^(z). What we can say safely is that the 
point density at z after taking an ensemble aver- 
age should be equal to the four-point function 

lim H 4 < e am e°" t '^e a ' l>iw) e " t>{z) >, 

w — >oo 

in the Liouville theory, where w is a point on the 
complex plane and a is the renormalization factor 
for the conformal mode. In a recent papcrM, 



an estimate for this four-point function has been 
given. This should be compared with Fig. 6. 
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